Spin-fermion model near the quantum critical point: one-loop renormalization group 
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We consider spin and electronic properties of itinerant electron systems, described by the spin- 
fermion model, near the antiferromagnetic critical point. We expand in the inverse number of hot 
spots in the Brillouin zone. A'' and present the results beyond previously studied N = oo limit. We 
found two new effects: (i) Fermi surface becomes nested at hot spots, and (ii) vertex corrections 
give rise to anomalous spin dynamics and change the dynamical critical exponent from z = 2 to 
z > 2. To first order in 1/iV we found z = 2N/{N - 2) which for a physical = 8 yields z ^ 2.67. 



PACS numbers: 74.20.Fg, 75.20Hr 

The problem of fermions interacting with critical an- 
tiferromagnetic spin fluctuations attracts a lot of atten- 
tion at the moment due to its relevance to both high 
temperature superconductors and heavy-fermion mate- 
rials The key interest of the current studies is to 
understand the system behavior near the quantum crit- 
ical point (QCP) where the magnetic correlation length 
diverges at T = |^] . Although in reality the QCP is al- 
most always masked by either superconductivity or pre- 
cursor effects to superconductivity, the vicinity of the 
QCP can be reached by varying external parameter such 
as pressure in heavy fermion compounds, or doping con- 
centration in cuprates. 

In this paper, we study the properties of the QCP with- 
out taking pairing fluctuations into account. We assume 
that the singularities associated with the closeness to the 
QCP extent up to energies which exceed typical energies 
associated with the pairing. This assumption is consis- 
tent with the recent calculations of the pairing instabil- 
ity temperature in cuprates From this perspective, 
the understanding of the properties of the QCP without 
pairing correlations is a necessary preliminary step for 
subsequent studies of the pairing problem. 

A detailed study of the antiferromagnetic QCP was 
performed by Hertz Q and later by Millis Q who chiefly 
focused on finite T properties near the QCP. They both 
argued that if the Fermi surface contains hot spots 
(points separated by antiferromagnetic momentum Q, 
see Fig. then spin excitations possess purely relax- 
ational dynamics with z = 2. They further argued that 
in d = 2, d + z = A, i.e., the critical theory is at marginal 
dimension, in which case one should expect that spin- 
spin interaction yields at maximum logarithmical cor- 
rections to the relaxational dynamics. Millis argued ||] 
that this is true provided that the effective Ginsburg- 
Landau functional for spins (obtained by integrating out 
the fermions) is an analytic function of the spin ordering 
field. This is a'priori unclear as the expansion coeffi- 
cients in the Ginsburg-Landau functional are made out 
of particle-hole bubbles and generally are sensitive to the 
closeness to quantum criticality due to feedback effect 



from near critical spin fiuctuations on the electronic sub- 
system. Millis however demonstrated that the quartic 
term in the Ginsburg-Landau functional is governed by 
high energy fermions and is free from singularities. 

In this communication, we, however, argue that the 
regular Ginsburg-Landau expansion is not possible in 2D 
by the reasons different from those displayed in [^,^. 
Specifically, we argue that the damping term in the spin 
propagator (assumed to be linear in uj in |^,|^) is by 
itself made out of a particle hole bubble, and, contrary 
to 0"* coefficient, is governed by low-energy fermions. We 
demonstrate that due to singular vertex corrections, the 
frequency dependence of the spin damping term at the 
QCP is actually w^^". In the one loop approximation, 
we find a w 0.25. 

Another issue which we study is the form of the renor- 
malized quasiparticle Fermi surface near the magnetic in- 
stability. In a mean-field SDW theory, the Fermi surface 
in a paramagnetic phase is not affected by the closeness 
to the QCP. Below the instability, the doubling of the 
unit cell induces a shadow Fermi surface at kp -\-Q, with 
the residue proportional to the deviation from critical- 
ity. This gives rise to the opening of the SDW gap near 
hot spots and eventually (for a perfect antiferromagnetic 
long range order) yields a Fermi surface in the form of 
small pockets around (7r/2,7r/2) and symmetry related 
points (see Fig. |l|a). Several groups argued |^ that this 
mean-field scenario is modified by fluctuations, and the 
Fermi surface evolution towards hole pockets begins al- 
ready within the paramagnetic phase. We show that the 
Fermi surface near hot spots does evolve as f ^ oo, but 
due to strong fermionic damping (not considered in [^), 
this evolution is a minor effect which at ^ = oo only gives 
rise to a nesting at the hot spots (see Fig. Qd). 

The point of departure for our analysis is the spin- 
fermion model which describes low-energy fermions inter- 
acting with their own collective spin degrees of freedom. 
The model is described by 
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FIG. 1. The Fermi surface with hot spots and the di- 
rections of Fermi velocities at hot spots separated by Q, and 
the evolution of the Fermi surface evolution for (a) mean-field 
(A'' = oo) SDW theory, and (b) finite A*'. In both cases, the 
doubling of the unit cell due to antiferromagnetic SDW order- 
ing introduces shadow Fermi surface and yields a gap opening 
near hot spots. At finite N , however, the Fermi surface at 
the quantum critical point becomes nested at hot spots due 
to vanishing of renormalized Vy. 
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Here cj^. ^ is the fermionic creation operator for an elec- 
tron with momentum k and spin projection a, cr^ are the 
Pauli matrices, and g measures the strength of the inter- 
action between fermions and their collective bosonic spin 
degrees of freedom. The latter are described by Sq and 
are characterized by a bare spin susceptibility which is 
obtained by integrating out high-energy fermions. 

This spin-fermion model can be viewed as the ap- 
propriate low-energy theory for Hubbard-type lattice 
fermion models provided that spin fluctuations are the 
only low-energy degrees of freedom. This model explains 
a number of measured features of cuprates both in the 
normal and the superconducting states [Q. Its applica- 
tion to heavy-fermion materials is more problematic as in 
these compounds conduction electrons and spins are in- 
dependent degrees of freedom, and the dynamics of spin 
fluctuations may be dominated by local Kondo physics 
rather than the interaction with fermions Q . 

The form of the bare susceptibility Xoil) is an input 
for the low-energy theory. We assume that xa (?) is non- 
singular and peaked at Q, i.e., Xo(q) = Xo/(C~^ + (q — 
Q)^), where ^ is the magnetic correlation length. In prin- 
ciple, xo can also contain a nonuniversal frequency de- 
pendent term in the form [uu/WY where W is of order 
of fermionic bandwidth. We, however, will see that for 
a Fermi surface with hot spots which we consider here, 
this term will be overshadowed by a universal cj^^" term 
produced by low-energy fermions. 

The earlier studies of the spin-fermion model have 
demonstrated that the perturbative expansion for both 
fermionic and bosonic self-energies holds in power of 
A = 3(7'^Xo/(47i'Ui?^~^) where vp is the Fermi velocity at 
a hot spot. This perturbation theory obviously does not 
converge when ^ ^ cx). As an alternative to a conven- 
tional perturbation theory, we suggested the expansion 



in inverse number of hot spots in the Brillouin zone N 
(= 8 in actual case) (|J^. Physically, large N implies 
that a spin fluctuation has many channels to decay into 
a particle-hole pair, which gives rise to a strong (~ N) 
spin damping rate. At the same time, a fermion near 
a hot spot can only scatter into a single hot spot sepa- 
rated by Q. Power counting arguments than show that 
a large damping rate appears in the denominators of the 
fermionic self-energy and vertex corrections and makes 
them small to the extent of The only exception 

from this rule is the fermionic self-energy due to a single 
spin fluctuation exchange, which contains a frequency de- 
pendent piece without prefactor due to an infrared 
singularity which has to be properly regularized |^ . 

The set of coupled equations for fermionic and bosonic 
self-energies at TV = oo has been solved in and we 
merely quote the result. Near hot spots, we have 



xo^V(i + (q-Q)'e'-*nn). 
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Here tk = Vxkx + Vyky, where k ^ k — khs, and v^, Vy, 
which we set to be positive, are the components of the 
Fermi velocity at a hot spot {vp = v'^+Vy). The fermionic 
self-energy Sfe(cL;) and the spin polarization operator Hfj 
are given by 
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and UJsf = (47r/iV) v^Vy/ig^xo^'^)- 

We see from Eq.(^ that for uj < uJsf, G{khs,^^) = 
Z/{uj+iuj\uj\/ (Aujsf)), i.e., as long as is finite, the system 
preserves the Fermi-liquid behavior at the lowest frequen- 
cies. The quasiparticle residue Z however depends on the 
interaction strength, Z = {1 + A)~^, and progressively 
goes down when the spin-fermion coupling increases. At 
larger frequencies to > uigf, the system crosses over to a 
region, which is in the basin of attraction of the quan- 
tum critical point, — oo. In this region, G^^{kF,oj) ~ 
3g {vxVyXo/T^NvlY/^ (?|w|)i/2sgn(w) §0. At the same 



time, spin propagator has a simple z = 2 relaxational dy- 
namics unperturbed by strong frequency dependence of 
the fermionic self-energy [ pd| . 

Our present goal is to go beyond N = oo limit and ana- 
lyze the role of 1/A^ corrections. The 1/N terms give rise 
to two new features: vertex corrections which renormal- 
ize both fermionic and bosonic self-energies, and static 
fermionic self-energy Efe. The corresponding diagrams 
are presented in Fig ^. The lowest-order 1 /N corrections 
have been calculated before Both vertex correc- 

tion and the static self-energy are logarithmical in ^: 
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straightforward, and yields 





FIG. 2. The one-loop RG diagrams for the fermionic 
self-energy and vertex renormalization. Solid lines are full 
fermionic propagators, wavy lines are full spin susceptibili- 
ties, and black triangles are full vertices. The lowest order 
diagrams are obtained by replacing full internal lines and ver- 
tices by their N = oo forms 

where ek+Q = -v^kx + Vyky, and Q{v) = 
(4/7r) arctan(t;2./t;j^) interpolates between Q = \ for Vj. = 
Vy, and Q = 2 for Vy 0. 

Besides, the corrections 

also contribute (l/-/V)wlog^ to G^^{lj), but this term 
is negligible compared to and we neglect it. 

We see from (^||) that the l/N corrections to the ver- 
tex and to the velocity of the excitations are almost de- 
coupled from each other: the velocity renormalization 
does not depend on the coupling strength at all, while the 
renormalization of the vertex depends on the ratio of ve- 
locities only through a non-singular Q{v). This is a direct 
consequence of the fact that the dynamical part of the 
spin propagator is obtained self-consistently within the 
model. Indeed, the overall factors in Ae^ and Ag/g are 
g'^iujsfS^'^) where oJafS,"^ comes from the dynamical part 
of the spin susceptibility. Since the fermionic damping 
is produced by the same spin-fermion interaction as the 
fermionic self-energy, uj^f scales as l/g^, and the coupling 
constant disappears from the r.h.s. of (^,^). 

The logarithmical dependence on ^ implies that 1/A^ 
expansion breaks down near the QCP, and one has to sum 
up the series of the logarithmical corrections. We will do 
this in a standard one-loop approximation by summing 
up the series in (l/iV)log^ but neglecting regular 
corrections to each term in the series. We verified that 
in this approximation, the cancellation of the coupling 
constant holds even when 5 is a running, scale dependent 
coupling. This in turn implies that one can separate the 
velocity renormalization from the renormalization of the 
vertex to all orders in 

Separating the corrections to and Vy and performing 
standard RG manipulations, we obtain a set of two RG 
equations for the running and 
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and Vy are the bare values of the 



where, we remind, v. 
velocities (the ones which appear in the Hamiltonian) . 

We see that v^ vanishes logarithmically at ^ — > 00. 
This implies that right at the QCP, the renormalized ve- 
locities at and fc/is + Q are antiparallel to each other, 
i.e. the Fermi surface becomes nested at hot spots (see 
Fig |l|b) . This nesting creates a "bottle neck effect" im- 
mediately below the criticality as the original and the 
shadow Fermi surfaces approach hot spots with equal 
derivatives (see Fig. |l|b). This obviously helps devel- 
oping a SDW gap at ky^g below the magnetic instability. 
However, above the transition, no SDW precursors ap- 
pear at T = 0. 

Another feature of the RG equations (||) is that they 
leave the product VxVy unchanged. This is a combination 
in which velocities appear in lOsj- The fact that v^Vy is 
not renormalized implies that, without vertex renormal- 
ization, uJsf^'^ remains finite at ^ = 00, i.e., spin fluctua- 
tions preserve a simple z = 2 relaxational dynamics. 

We now consider vertex renormalization. Using again 
the fact that g'^uJsf does not depend on the running 
coupling constant, one can straightforwardly extent the 
second-order result for the vertex renormalization, Eqn 
(0), to the one-loop RG equation 
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where g^ is a running coupling constant, and Q{v) is the 
same as in (jj) but contain renormalized velocities v^ and 
Vy. At the QCP, the dependence on ^ obviously trans- 
forms into the dependence on frequency (L = log^ — > 
(1/2) log jo^o/ti^l, where ujq is the upper cutoff). Using the 
fact that for ^ — > 00, v^/v^- w Ntt/24L and expanding 
Q{v) near = 0, we find Q{v) « 2(1 - (2/tt)v^/v^) = 
2 — N/3L. Substituting this result into (||) and solving 
the differential equation we obtain {lu = uj/luq) 
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We see that at the QCP, running coupling constant di- 
verges as a; ^ roughly as |a;|~"^/^. Substituting this 
result into the spin polarization operator and using the 
fact that ojsf oc ((?^)~^ we find that at the QCP, 
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where L = log^. The solution of these equations is 



This result implies that vertex corrections change the dy- 
namical exponent z from its mean-field value z = 2 to 
z = 2N/{N - 2). For = 8, this yields z « 2.67 and 
x{Q,ijj) oc where a = 0.75. 



3 



Singular vertex corrections also renormalize the 
fermionic self-energy as oc g^^J\uj\/vF- Using the 

results for and vp ^ Vx ^e. obtain at criticality 



logw| 



(11) 



Eqs. (0), ([lo| ) and ( p]) are the central results of the pa- 
per. We see that the singular corrections to the Fermi 
velocity cause nesting but do not affect the spin dynam- 
ics. The corrections to the vertex on the other hand do 
not affect velocities, but change the dynamical critical 
exponent for spin fluctuations. 

We now briefly discuss the form of the susceptibility 
at finite T. Previous studies have demonstrated [^|J|] 
that the scattering of a given spin fluctuation by classi- 
cal, thermal spin fluctuations yields, up to logarithmical 



prefactors, ^ 



uT, where u is the coefficient in the 



(j)^ term in the Ginsburg-Landau potential. This implies 
that at the QCP, x((9,cj) (xT -i\uj\. 

We, however, argue that the linear in T and the linear 
in uj terms have completely different origin: the linear in 
Lo term comes from low-energies and is universal, while 
the linear in T term comes from high energies and is 
model dependent. This can be understood by analyzing 
the particle-hole bubble at finite T. We found that as 
long as one restricts with the linear expansion near the 
Fermi surface. Ho preserves exactly the same form as at 
T = 0, to all orders in the perturbation theory. The tem- 
perature dependence of 11 appears only due to a nonzero 
curvature of the electronic dispersion and is obviously 
sensitive to the details of the dispersion at energies com- 
parable to the bandwidth. Similarly, the derivation of 
the Landau-Ginsburg potential from shows Q that 
u vanishes for linearized e^, and is finite only due to a 
nonzero curvature of the fermionic dispersion. 

The different origins of T and uj dependences in x{Q^ ^) 
imply that the anomalous frequency dependence of 
xiQ,^) is not accompanied by the anomalous tempera- 
ture dependence of x('5; 0) simply because for high en- 
ergy fermions, vertex corrections are non-singular. This 
result implies, in particular, that our theory does not 
explain anomalous spin dynamics observed in heavy 
fermion p3| despite the similarity in the exponent for 
the frequency dependence of IIq, because the experi- 
mental data imply the existence of the O/T scaling in 
CeCu(,-xAux ■ More likely, the explanation should in- 
volve the local Kondo physics |8| . 

Finally, we consider how anomalous vertex corrections 
affect the superconducting problem. We and Finkel'stein 
argued recently § that at ^ = oo, the kernel K{uj,Vt) of 
the Eliashberg-type gap equation for the d— wave anoma- 
lous vertex F{Vl) = (7rr/2) Y,u> ^i^^ ^)F{uj) behaves as 
K{uj,fl) oc gV(t,2,s2(a;)nn-^)i/2 At iV = oo, this yields 
(including prefactor) — ~ u})\^^^'^ . Al- 

though this kernel is qualitatively different from the one 
in the BCS theory because it depends on both frequen- 



cies, it still scales as inverse frequency due to an inter- 
play between a non-Fermi liquid form of the fermionic 
self-energy and the absence of the gap in the spin sus- 
ceptibility which mediates pairing. We demonstrated in 
1^ that this inverse frequency dependence gives rise to a 
finite pairing instability temperature even when ^ = oo. 

To check how the kernel is affected by vertex correc- 
tions, we substitute the results for g^, vp, and Iln 
into K{Lu,n). We find after simple manipulations that 
despite singular vertex corrections, the kernel in the gap 
equation still scales inversely proportional to frequency. 
A simple extension of the analysis in |^ then shows that 
the system still possesses a pairing instability at ^ = oo 
at a temperature which differs from that without vertex 
renormalization only by 1/iV corrections. 

To summarize, in this paper we considered the prop- 
erties of the antiferromagnetic quantum critical point for 
itinerant electrons by expanding in the inverse number of 
hot spots in the Brillouin zone TV = 8. We went beyond a 
self-consistent N — oo theory and found two new effects: 
(i) Fermi surface becomes nested at hot spots which is 
a weak SDW precursor effect, and (ii) vertex corrections 
account for anomalous spin dynamics and change the dy- 
namical critical exponent from z — 2 to z > 2. To first 
order in 1/N we found z = 2N/{N - 2) Ki 2.67. We 
argued that anomalous frequency dependence is not ac- 
companied by anomalous T dependence. 
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